Spin polarized states in neutron matter at a strong magnetic field 
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Spin polarized states in neutron matter at strong magnetic fields up to 10^* G are considered 
in the model with the Skyrme effective interaction. By analyzing the self-consistent equations at 
zero temperature, it is shown that a thermodynamically stable branch of solutions for the spin 
polarization parameter as a function of density corresponds to the negative spin polarization when 
the majority of neutron spins are oriented opposite to the direction of the magnetic field. Besides, 
beginning from some threshold density dependent on the magnetic field strength the self-consistent 
equations have also two other branches of solutions for the spin polarization parameter with the 
positive spin polarization. The free energy corresponding to one of these branches turns out to 
be very close to that of the thermodynamically preferable branch. As a consequence, at a strong 
magnetic field, the state with the positive spin polarization can be realized as a metastable state at 
the high density region in neutron matter which under decreasing density at some threshold density 
changes into a thermodynamically stable state with the negative spin polarization. 
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I. INTRODUCTION 

Neutron stars observed in nature are magnetized ob- 
jects with the magnetic field strength at the surface in 
the range of 10^-10^^ G For a special class of neutron 
stars such as soft gamma-ray repeaters and anomalous 
X-ray pulsars, the field strength can be much larger and 
is estimated to be about lO^^'-lO^^ G 0. These strongly 
magnetized objects are called magnetars 0] and comprise 
about 10% of the whole population of neutron stars [I]- 
However, in the interior of a magnetar the magnetic field 
strength may be even larger, reaching values of about 
10^* G [1, Q. The possibility of existence of such ultra- 
strong magnetic fields is not yet excluded, because what 
we can learn from the magnetar observations by their pe- 
riods and spin-down rates, or by hydrogen spectral lines 
is only their surface fields. There is still no general con- 
sensus regarding the mechanism to generate such strong 
magnetic fields of magnetars, although different scenarios 
were suggested such as, e.g., a turbulent dynamo ampli- 
fication mechanism in a neutron star with the rapidly ro- 
tating core at first moments after it goes supernova [2] , or 
the possibility of spontaneous spin ordering in the dense 
quark core of a neutron star Q . 

Under such circumstances, the issue of interest is the 
behavior of neutron star matter in a strong magnetic 
field d, H, H, 0. In the recent study 0], neutron star 
matter was approximated by pure neutron matter in 
the model considerations with the effective Skyrme and 



Gogny forces. It has been shown that the behavior of the 
spin polarization of neutron matter in the high density 
region at a strong magnetic field crucially depends on 
whether neutron matter develops a spontaneous spin po- 
larization (in the absence of a magnetic field) at several 
times nuclear matter saturation density as is usual for the 
Skyrme forces, or the appearance of a spontaneous polar- 
ization is not allowed at the relevant densities (or delayed 
to much higher densities), as in the case with the Gogny 
DIP force. In the former ferromagnetic transi- 

tion to a totally spin polarized state occurs while in the 
latter case a ferromagnetic transition is excluded at all 
relevant densities and the spin polarization remains quite 
low even in the high density region. Note that the issue of 
spontaneous appearance of spin polarized states in neu- 
tron and nuclear matter is a controversial one. On the 
one hand, the models with the Skyrme effective nucleon- 
nucleon (NN) interaction predict the occurrence of spon- 
taneous spin instability in nuclear matter at densities in 
the range from qq to 4pn for different parametrizations of 
the NN potential [lo|-|23| (go = 0.16 fm"^ is the nuclear 
saturation density). For the Gogny effective interaction, 
a ferromagnetic transition in neutron matter occurs at 
densities larger than Iqq for the DIP parametrization 
and is not allowed for Dl, DIS parametrizations [23j . 
However, for the DIS Gogny force an antiferromagnetic 
phase transition happ ens in symmetric nuclear matter at 
the density S.Sgo [2J|- On the other hand, for the models 
with the realistic NN interaction, no sign of spontaneous 
spin instability has been found so far at any isospin asym- 
metry up to densities well above qq [2^-[3i[. 
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Here we study thermodynamic properties of spin po- 
larized neutron matter at a strong magnetic field in the 
model with the Skyrme effective forces. As a framework 
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for consideration, we choose a Fermi liquid approach for 
the description of nuclear matter [H, [H, [s^l • Proceeding 
from the minimum principle for the thermodynamic po- 
tential, we get the self-consistent equations for the spin 
order parameter and chemical potential of neutrons. In 
the absence of a magnetic field, the self-consistent equa- 
tions have two degenerate branches of solutions for the 
spin polarization parameter corresponding to the case, 
when the majority of neutron spins are oriented along, 
or opposite to the spin quantization axis (positive and 
negative spin polarization, respectively). In the presence 
of a magnetic field, these branches are modified differ- 
ently. A thermodynamically stable branch corresponds 
to the state with the majority of neutron spins aligned 
opposite to the magnetic field. At a strong magnetic 
filed, the branch corresponding to the positive spin po- 
larization splits into two branches with the positive spin 
polarization as well. The last solutions were missed in the 
study of Ref. 0. We perform a thermodynamic analysis 
based on the comparison of the respective free energies 
and arrive at the conclusion about the possibility of the 
formation of metastable states in neutron matter with 
the majority of neutron spins directed along the strong 
magnetic field. The appearance of such metastable states 
can be possible due to the strong spin-dependent medium 
correlations in neutron matter with the Skyrme forces at 
high densities. 

Note that we consider thermodynamic properties of 
spin polarized states in neutron matter at a strong mag- 
netic field up to the high density region relevant for astro- 
physics. Nevertheless, we take into account the nucleon 
degrees of freedom only, although other degrees of free- 
dom, such as pions, hyperons, kaons, or quarks could be 
important at such high densities. 



II. BASIC EQUATIONS 



Y4 = — /.jq/T being the Lagrange multipliers, fio being the 
chemical potential of neutrons, and T the temperature. 

Given the possibility for alignment of neutron spins 
along or opposite to the magnetic field H, the normal 
distribution function of neutrons and single particle en- 
ergy can be expanded in the Pauli matrices ai in spin 
space 



/(P) = /o(p)cro + /3(p)cr3, 

^(P) = eo(p)cro + e3(p)a•3■ 



(3) 



Using Eqs. ([2]) and ([3]), one can express evidently the 
distribution functions /o, fs in terms of the quantities e: 



fo = ^{n(uj+) +n{uj-)}, 
h = -n{uj^)}. 

Here n{Lo) — {exp(ybw) 4- 1}~^ and 

^± = Co ±6, 
Co = £0 — A*o, C3 = £3- 



(4) 



(5) 



As follows from the structure of the distribution func- 
tions /, the quantity being the exponent in the Fermi 
distribution function n, plays the role of the quasiparticle 
spectrum. The spectrum is twofold split due to the spin 
dependence of the single particle energy ^(p) in Eq. ([3]). 
The branches lj± correspond to neutrons with spin up 
and spin down. 

The distribution functions / should satisfy the norma- 
lization conditions 



^E/o(p) 



V 



^fsiv) = Q}- Q[= Ap. 



(6) 
(7) 



The normal (nonsuperfluid) states of neutron matter 
are described by the normal distribution function of neu- 
trons /kik2 — Tr Qo^^ati-^^, where k = (p, cr), p is momen- 
tum, a is the projection of spin on the third axis, and g 
is the density matrix of the system [2l| . Further it 
will be assumed that the third axis is directed along the 
external magnetic field H. The energy of the system is 
specified as a functional of the distribution function /, 
E = E{f), and determines the single particle energy 



dEjf) 



(1) 



The self-consistent matrix equation for determining the 
distribution function / follows from the minimum condi- 
tion of the thermodynamic potential [s^ . [s^ and is 



/ = {exp(yo£ + n) + 1}"' = {exp(yoe) + 1}" 



(2) 



Here the quantities e and I4 are matrices in the space 
of K variables, with Y^i^-^^^ = Y4Sk,-i^k2, Yq = 1/T, and 



Here g = + Qi is the total density of neutron matter, 
gi and gi are the neutron number densities with spin 
up and spin down, respectively. The quantity Ag may 
be regarded as the neutron spin order parameter. It de- 
termines the magnetization of the system M = fijiAg, 
fi„ being the neutron magnetic moment. The magne- 
tization may contribute to the internal magnetic field 
B = H + AttM . However, we will assume, analogously to 
Refs. 0, that the contribution of the magnetization 
to the magnetic field B remains small for all relevant 
densities and magnetic field strengths, and, hence. 



B ^ H. 



(8) 



This assumption holds true due to the tiny value of the 
neutron magnetic moment /i„ = — 1.9130427(5)/ijv ~ 
-6.031 • 10-18 MeV/G [H {fiN being the nuclear magne- 
ton) and is confirmed numerically by finding solutions of 
the self-consistent equations in two approximations, cor- 
responding to preserving and neglecting the contribution 
of the magnetization. 
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In order to get the self-consistent equations for the 
components of the single particle energy, one has to set 
the energy functional of the system. In view of the ap- 
proximation dS]) , it reads [2l|, [s^l 



Eo{f,H) = 2^eo(p)/o(p) - 2^„i/^/3(p), 



(9) 



Eintif) 



E 



field 



E{^~0(P)/0(P)+=~3(P)/3(P)}, 

P 
IP 



V, 



where 



eoip) = ^^C/o"(k)/o(q), k= (10) 
q 

e~3(p) = ^Et^r(k)/3(q)- 

q 

2 

Here eo(p) = is the free single particle spectrum, 
Too is the bare mass of a neutron, ?/o''(k), ?7" (k) are the 
normal Fermi liquid (FL) amplitudes, and £0,23 are the 
FL corrections to the free single particle spectrum. Note 
that in this study we will not be interested in the total 
energy density and pressure in the interior of a neutron 
star. By this reason, the field contribution £'/ie;d, being 
the energy of the magnetic field in the absence of matter, 
can be omitted. Using Eqs. ([T|) and 0, we get the self- 
consistent equations in the form 



Co(p) = £o(p) + £q{p) - /^o, 6(P) 



(11) 



To obtain numerical results, we utilize the effective 
Skyrme interaction. The amplitude of NN interaction 
for the Skyrme effective forces reads [s^ 



?}(p,q) = to{l + XoPa) 

1 



1 



2h? 



<2 



(12) 

(1 4- X2Pa)pCl, 



where Po- = (1 + (Ti(T2)/2 is the spin exchange opera- 
tor, ti,Xi and /3 are some phenomenological parameters 
specifying a given parametrization of the Skyrme inter- 
action. In Eq. p2)) . the spin-orbit term irrelevant for a 
uniform matter was omitted. The normal FL amplitudes 
can be expressed in terms of the Skyrme force parame- 
ters [33,0: 



t^o(k) 



t^r(k) 



2to{l~xo) + jg''{l-X3) 
2 

-2to(l 
2 



Xq 



Xi) + 3t2{l + X2W 



)-^/(l-^3) 



Ml + X2) - hii - xi)]k 



2 



(13) 

(14) 
6„k2. 



Further we do not take into account the effective tensor 
forces, which lead to coupling of the momentum and spin 
degrees of freedom i37., and, correspondingly, to 

anisotropy in the momentum dependence of FL ampli- 
tudes with respect to the spin quantization axis. Then 



^0 



2mr. 



6 = -fJ-nH + {an + b,i—) ^ 
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(15) 
(16) 



where the effective neutron mass m„ is defined by the 
formula 



9 



2to„ 



2mo 8 



+ ^[<i(l-xi)+ 3*2(1 + ^2)], (17) 



and the renormalized chemical potential fi should be de- 
termined from Eq. ©. The quantity (q^)3 in Eq. p6p is 
the second order moment of the distribution function f^: 



q')3 = -^Eq'/3(q). 



(18) 



In view of Eqs. (|T5|) . p6)) . the branches uj± = of the 
quasiparticle spectrum in Eq. ([5]) read 

M + a(-^j„i/ + + — (q^)3), (19) 



2to„ 



where to^ is the effective mass of a neutron with spin up 
(cr = -|-1) and spin down (cr = —1) 



2ma 



2mo 



+ ^t2{l+X2) 



(20) 



-^hil - a;i) + ^2(1 + X2)], £»+(-) = (?T(i)- 



Note that for totally spin polarized neutron matter 
Too . , grriQ 



— = 1 -f- -r^t2(l + X2), 
TO* 



(21) 



where to* is the effective neutron mass in the fully po- 
larized state. Since usually for Skyrme parametrizations 
^2 < 0, we have the constraint X2 < —1, which guaran- 
tees the stability of totally polarized neutron matter at 
high densities. 

It follows from Eq. (fTO|) that the effective chemical po- 
tential /iCT for neutrons with spin- up (cr = 1) and spin- 
down ((T = — 1) can be determined as 



16 ' 



a(M„iJ- = -^(q^)3) 



(22) 



Thus, with account of expressions ^ for the distri- 
bution functions /, we obtain the self-consistent equa- 
tions (O , ([7]) , and (fT8|) for the effective chemical potential 
/i, spin order parameter A^?, and second order moment 

(q')3. 
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III. SOLUTIONS OF SELF-CONSISTENT 
EQUATIONS AT T = 0. THERMODYNAMIC 
STABILITY 

Here we directly solve the self-consistent equations at 
zero temperature and present the neutron spin order 
parameter as a function of density and magnetic field 
strength. In solving numerically the self-consistent equa- 
tions, we utilize SLy4 and SLy7 Skyrme forces [i^l, which 
were constrained originally to reproduce the results of 
microscopic neutron matter calculations (pressure versus 
density curve) . Note that the density dependence of the 
nuclear symmetry energy, calculated with these Skyrme 
interactions, gives the neutron star models in a broad 
agreement with the obscrvablcs such as the minimum 
rotation period, gravitational mass-radius relation, the 
binding energy, released in supernova collapse, etc. [4]| . 
Besides, these Skyrme parametrizations satisfy the con- 
straint X2 < —1, obtained from Eq. (PT|) . 

We consider magnetic fields up to the values allowed 
by the scalar virial theorem. For a neutron star with 
the mass M and radius R, equating the magnetic field 
energy Eh ^ {inR^ /3){H^ /Sn) with the gravitational 
binding energy Eq ^ GJVP/R, one gets the estimate 
Umax ~ -^(SG)^/^. For a typical neutron star with M — 
1.5Mq and R = 1Q~^Rq, this yields for the maximum 
value of the magnetic field strength Hmax ~ 10^* G. This 
magnitude can be expected in the interior of a magnctar 
while recent observations report the surface values up 
to H ^ 10^^ G, as inferred from the hydrogen spectral 
lines m. 

In order to characterize spin ordering in neutron mat- 
ter, it is convenient to introduce a neutron spin polariza- 
tion parameter 



n 



Q 



Ag 

Q 



(23) 



Fig. [T] shows the dependence of the neutron spin po- 
larization parameter from density, normalized to the nu- 
clear saturation density go, at zero temperature in the 
absence of the magnetic field. The spontaneous polar- 
ization develops at g = S.TOgo for the SLy4 interaction 
(go ~ 0.16 fm^"^) and at g ~ 3.59£)o for the SLy7 in- 
teraction {go — 0.158 fm^'^), that reflects the instabil- 
ity of neutron matter with the Skyrme interaction at 
such densities against spin fluctuations. Since the self- 
consistent equations at iJ = are invariant with respect 
to the global flip of neutron spins, we have two branches 
of solutions for the spin polarization parameter, nj(p) 
(upper) and 11,^ (^?) (lower) which differ only by sign, 
Il+ig)^-U^{g). 

Fig. [5] shows the neutron spin polarization parameter 
as a function of density for a set of fixed values of the 
magnetic field. The branches of spontaneous polariza- 
tion are modified by the magnetic field differently, since 
the self-consistent equations at i7 ^ lose the invari- 
ance with respect to the global flip of the spins. At 
nonvanishing H, the lower branch Ili{g), corresponding 
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FIG. 1: (Color online) Neutron spin polarization parameter as 
a function of density at vanishing temperature and magnetic 
field. 



to the negative spin polarization, extends down to the 
very low densities. There are three characteristic den- 
sity domains for this branch. At low densities g < Q.hgo-, 
the absolute value of the spin polarization parameter in- 
creases with decreasing density. At intermediate densi- 
ties 0.5go 1^ SgiQ, there is a plateau in the ni(gi) de- 
pendence, whose characteristic value depends on iJ, e.g.. 
Hi « -0.08 at 7? = IQis G. At densities g > 3go, the 
magnitude of the spin polarization parameter increases 
with density, and neutrons become totally polarized at 

g « e^o- 

Note that the results in the low-density domain should 
be considered as a first approximation to the real complex 
picture, since, as discussed in detail in Ref. @, the low 
density neutron-rich matter in /3-equilibrium possesses 
a frustrated state, "nuclear pasta", arising as a result 
of competition of Coulomb long-range interactions and 
nuclear short-range forces. In our case, where a pure 
neutron matter is considered, there is no mechanical in- 
stability due to the absence of the Coulomb interaction. 
However, the possibility of appearance of low-density nu- 
clear magnetic pasta and its impact on the neutrino opac- 
ities in the protoneutron star early cooling stage should 
be explored in a more detailed analysis. 

Let us consider the modification of the upper branch 
of spontaneous polarization Hj(gi) at nonvanishing mag- 
netic field. It is seen from Fig.[2]that now beginning from 
some threshold density the self-consistent equations at a 
given density have two positive solutions for the spin po- 
larization parameter (apart from one negative solution). 
These solutions belong to two branches, Il2{g) and 113(g), 
characterized by different dependence from density. For 
the branch H2(g), the spin polarization parameter de- 
creases with density and tends to zero value while for the 
branch Tisig) it increases with density and is saturated. 
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FIG. 2: (Color online) Neutron spin polarization parameter 
as a function of density at T = and different magnetic field 
strengths for (a) SLy4 interaction and (b) SLy7 interaction. 
The branches of spontaneous polarization IIq , IIq' are shown 
by solid curves. 
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FIG. 3: (Color online) Energy per neutron as a function of 
density at T = for different branches Ili{g)-'n3{g) of solu- 
tions of the self-consistent equations at H = 10^* G for (a) 
SLy4 and (b) SLy7 interactions, including a spontaneously 
polarized state. 



These branches appear step-wise at the same threshold 
density pth dependent on the magnetic field and being 
larger than the critical density of spontaneous spin in- 
stability in neutron matter. For example, for SLy7 inter- 
action, pth ~ 3.80 go a.t H = 5 - 10^'^ G, and gth ~ 3.92 go 
dX H = 10^® G. The magnetic field, due to the negative 
value of the neutron magnetic moment, tends to orient 
the neutron spins opposite to the magnetic field direc- 
tion. As a result, the spin polarization parameter for the 
branches 112 (gi), ^^{g) with the positive spin polariza- 
tion is smaller than that for the branch of spontaneous 
polarization 11^, and, vice versa, the magnitude of the 
spin polarization parameter for the branch T\i{g) with the 
negative spin polarization is larger than the correspond- 
ing value for the branch of spontaneous polarization IIq" . 
Note that the impact of even such strong magnetic field 
as H ^ 10^^ G is small: The spin polarization parameter 
for all three branches Iii{g)-Ii-i{Q) is either close to zero, 
or close to its value in the state with spontaneous polar- 



ization, which is governed by the spin-dependent medium 
correlations. 

Thus, at densities larger than gth, we have three 
branches of solutions: one of them, T\i{g), with the neg- 
ative spin polarization and two others, 112 (p) and Ii^[g), 
with the positive polarization. In order to clarify, which 
branch is thermodynamically preferable, we should com- 
pare the corresponding free energies. Fig. [3] shows the 
energy per neutron as a function of density at T = and 
H = 10^® G for these three branches, compared with the 
energy per neutron for a spontaneously polarized state 
[the branches 11^ (p)]. It is seen that the state with the 
majority of neutron spins oriented opposite to the direc- 
tion of the magnetic field [the branch ni(p)] has a lowest 
energy. This result is intuitively clear, since magnetic 
field tends to direct the neutron spins opposite to H, 
as mentioned earlier. However, the state, described by 
the branch n3(g) with the positive spin polarization, has 
the energy very close to that of the thermodynamically 
stable state. This means that despite the presence of a 
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FIG. 4: (Color online) Spin polarization parameter as a 
function of the magnetic field strength at T = for differ- 
ent branches ni(_H')-n3(_ff) of solutions of the self-consistent 
equations at £i = Aqo and for the branch ni(_ff) at £i = 2qo 
for (a) SLy4 interaction and (b) SLy7 interaction. 



strong magnetic field H ~ 10^® G, the state with the 
majority of neutron spins directed along the magnetic 
field can be realized as a metastable state in the dense 
core of a neutron star in the model consideration with 
the Skyrmc effective interaction. In this scenario, since 
such states exist only at densities g ^ pth , under decreas- 
ing density (going from the interior to the outer regions 
of a magnetar) a metastable state with the positive spin 
polarization at the threshold density pth changes into a 
thermodynamically stable state with the negative spin 
polarization. 

At this point, note some important differences be- 
tween the results in our study and those obtained in 
Ref. First, in the study [9| of neutron matter at a 
strong magnetic field only one branch of solutions for 
the spin polarization parameter was found in the model 
with the Skyrme interaction (for the same SLy4 and 
SLyT parametrizations). However, in fact, we have seen 
that the degenerate branches of spontaneous polarization 
(at zero magnetic field) with the positive and negative 



FIG. 5: (Color online) Same as in Fig. |4] but for the energy 
per neutron. 



spin polarization are modified differently by the mag- 
netic field, and, as a result, in the Skyrme model, in 
general, there are three different branches of solutions 
of the self-consistent equations at nonvanishing magnetic 
field. Besides, the only branch considered in Ref. [3 and 
corresponding to our thermodynamically stable branch 
Hi, is characterized by the positive spin polarization, 
contrary to our result with Hi < 0. This disagreement 
is explained by the incorrect sign before the term with 
the magnetic field in the equation for the quasiparticlc 
spectrum in Ref. Q (analogous to Eq. ([TO]) in our case). 
Clearly, in the equilibrium configuration the majority of 
neutron spins are aligned opposite to the magnetic field. 

Fig.[3]shows the spin polarization parameter as a func- 
tion of the magnetic field strength at zero temperature 
for different branches Ili{H)-Il3{H) of solutions of the 
self-consistent equations at ^3 = Ago compared with that 
for the branch Ili{H) a,t g = 2go. It is seen that up to 
the field strengths H = 10^^ G, the infiuence of the mag- 
netic field is rather marginal. For the branches Ili{H) 
and 112(H), the magnitude of the spin polarization pa- 
rameter increases with the field strength while for the 
113(H) it decreases. Interestingly, as is clearly seen from 
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the top panel for the SLy4 interaction, at the given den- 
sity, there exists some maximum magnetic field strength 
-ffm at which the branches 112 a-nd II^ converge and do 
not continue at H > H„i. 

Fig. [5] shows the energy of neutron matter per particle 
as a function of the magnetic field strength at T = un- 
der the same assumptions as in Fig. ID It is seen that the 
state with the negative spin polarization [branch Ili(H)] 
becomes more preferable with increasing the magnetic 
field although the total effect of changing the magnetic 
field strength by two orders of magnitude on the en- 
ergy corresponding to all three branches ni(_ff)-n3(_ff) 
remains small. It is also seen that the increase of the 
density by a factor of two leads to the increase in the en- 
ergy per neutron roughly by a factor of three reflecting 
the fact that the medium correlations play more impor- 
tant role in building the energetics of the system than 
the impact of a strong magnetic field. 

IV. CONCLUSIONS 

We have considered spin polarized states in neu- 
tron matter at a strong magnetic field in the model 
with the Skyrme effective NN interaction (SLy4, SLy7 
parametrizations). The self-consistent equations for the 
spin polarization parameter and chemical potential of 
neutrons have been obtained and analyzed at zero tem- 
perature. It has been shown that the thermodynamically 
stable branch of solutions for the spin polarization pa- 
rameter as a function of density corresponds to the case 
when the majority of neutron spins are oriented oppo- 
site to the direction of the magnetic field (negative spin 
polarization). This branch extends from the very low 
densities to the high density region where the spin po- 
larization parameter is saturated, and, respectively, neu- 
trons become totally spin polarized. Besides, beginning 
from some threshold density gth being dependent on the 
magnetic field strength the self-consistent equations have 
also two other branches (upper and lower) of solutions for 
the spin polarization parameter corresponding to the case 
when the majority of neutron spins are oriented along the 
magnetic field (positive spin polarization). For example, 
for SLy7 interaction, gth ~ 3.80 go a.t H = 5 ■ 10^^ G, 
and Qth ~ 3.92 qq at H = 10^® G. The spin polariza- 
tion parameter along the upper branch increases with 
density and is saturated, while along the lower branch it 
decreases and vanishes. The free energy corresponding to 
the upper branch turns out to be very close to the free en- 
ergy corresponding to the thermodynamically preferable 
branch with the negative spin polarization. As a con- 
sequence, at a strong magnetic field, the state with the 
positive spin polarization can be realized as a metastable 
state at the high density region in neutron matter which 



under decreasing density (going from the interior to the 
outer regions of a magnetar) at the threshold density gth 
changes into a thermodynamically stable state with the 
negative spin polarization. 

In this study, we have considered the zero tempera- 
ture case, but as was shown in Ref. 0], the influence of 
finite temperatures on spin polarization remains moder- 
ate in the Skyrme model, at least, up to the temperatures 
relevant for protoneutron stars, and, hence, one can ex- 
pect that the considered scenario will be preserved at 
finite temperatures as well. The possible existence of a 
metastable state with positive spin polarization will af- 
fect the neutrino opacities of a neutron star matter in a 
strong magnetic field, and, hence, will lead to the change 
of cooling rates of a neutron star compared to cooling 
rates in the scenario with the majority of neutron spins 
oriented opposite to the magnetic field [ist . 

The calculations of the neutron spin polarization pa- 
rameter and energy per neutron show that the infiuence 
of the magnetic field remains small at the field strengths 
up to 10^^ G. Note that in Ref. [H the consideration 
also has been done for the Gogny effective NN interac- 
tion (DIS, DIP parametrizations) up to densities 4po- 
Since for the DIS parametrization there is no sponta- 
neous ferromagnetic transition in neutron matter for all 
relevant densities, and for the DIP parametrization this 
transition occurs at the density larger than 7 go [23j . no 
sign of a ferromagnetic transition at a strong magnetic 
field was found in Ref. up to densities Ago for these 
Gogny forces. According to our consideration, one can 
expect that the metastable states with the positive spin 
polarization in neutron matter at a strong magnetic field 
could appear at densities larger than 7go for the DIP 
parametrization while the scenario with the only branch 
of solutions corresponding to the negative spin polariza- 
tion would be realized for the DIS force. 

It is worthy to note also that in the given research a 
neutron star matter was approximated by pure neutron 
matter. This approximation allows one to get the qual- 
itative description of the spin polarization phenomena 
and should be considered as a first step towards a more 
realistic description of neutron stars taking into account 
a finite fraction of protons with the charge neutrality and 
beta equilibrium conditions. In particular, some admix- 
ture of protons can affect the onset densities of enhanced 
polarization in a neutron star matter with the Skyrme 
interaction. 
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